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Abstract 

Thermoelectric power of a material, typicahy governed by its band structure and carrier density, 
can be varied by chemical doping that is often restricted by solubility of the dopant. Materials 
showing large thermoelectric power are useful for many industrial applications, such as the heat- 
to-electricity conversion and the thermoelectric cooling device. Here we show a full electric field 
tuning of thermoelectric power in a dual-gated bilayer graphene device resulting from the opening 
of a band-gap by applying a perpendicular electric field on bilayer graphene. We uncover a large 
enhancement in thermoelectric power at low temperature, which may open up a new possibility in 
low temperature thermoelectric application using graphene-based device. 

PACS numbers: 65.80.Ck, 73.63.Bd, 72.15.Jf, 73.63.-b 
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Bilayer graphene (BLG) comprises two monolayers of graphene stacked as in bulk graphite 
exhibiting unusual massive chiral fermionic excitations [l|-l3|. Its conduction band and va- 
lence band touch at charge neutral point (CNP) resembling a zero-gap semiconductor with 
an effective mass m* ^ 0.054me. Many interests reside on the band gap engineering in BLG 
by breaking its inversion symmetry, which is of particular importance for semiconductor 
device application. Ohta et.al. jj] demonstrated the band gap opening by chemical doping 
of potassium onto the upper layer of a BLG placed on a silicon carbide substrate using 
angle-resolved photoemission spectroscopy. Soon later, Castro et.al. showed the band 
gap tuning by electric field effect in a BLG device on SiO2(300 nm)/Si substrate with NH3 
adsorption on the upper layer to further enhance the layer asymmetry. More recently, a 
dual-gated geometry for band gap engineering was realized in BLG devices 6|-[8|, where a 
BLG was sandwiched by top gate and bottom gate providing a full electric field control over 
its band structure. However, careful treatment for the top gate dielectric turns out to be 
crucial in order to avoid the downfall of mobility in BLG 9|, Il0| • 

In this report, we utilize high quality micro-crystals of hexagonal boron nitride (hBN) as 
;op gate dielectric, which was shown to impose less trapped charges and strain on graphene 



The dual-gated BLG devices we fabricated [12] show carrier mobility (~ 2,000-3,000 



cm^/V-sec) comparable to the ones without top-gate/hBN. We demonstrate a full electric- 
field tuning of thermoelectric power (TEP) in bilayer graphene devices, which has been 
predicted theoretically [3] and originates from the opening of a band gap via the application 
of an out-of-plane electrical displacement field D. We uncover an enhancement in TEP due 
to D, which grows larger at lower temperatures attaining a value comparable to or exceeding 
those of Bi2Te3-based alloys and sodium cobaltates NaxCo02 at similar temperature. Our 
result reveals the potential thermoelectric application using graphene-based device. 

The device geometry, which comprises a heater, two local thermometers (Rtl and Rt2), 
three voltages leads and a local t op g ate as shown in Fig. [I](a), enables the 4-probe mea- 



surements of resistance and TEP [12, 



14 



15| . The thickness of the hBN micro-crystal for 



this particular device is ~ 40 nm determined by an atomic force microscope. We note that 
the top-gate effective region, shown as shaded area in Fig. [T](a), only partially covers the 
BLG area between voltage leads. Therefore, additional geometric factor needs to be taken 
into account in order to extract the sheet resistance and TEP for BLG under the top-gate 
effective region |l2|. For the following discussion, we used the suffix e(ne) referring to the 
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top-gate (non)effective region and suffix □ for the sheet resistance. 

Figure lH^b) shows a contour plot for the measured resistance R at 200 K for the device 
shown in Fig. dJ^a). The bottom-gate vohage Vbg is swept up to ± 70 V, while the top-gate 
voltage Vtg is kept at certain value from -10 V to lOV. The high R states occur in regions 
with large Vbg and Vtg at opposite polarity defined in Fig. [U^a). When plotting (Vtg, Vbg) for 
the peak position of R (Fig. [U^d)), it is nearly temperature independent and shows excellent 
linearity giving a slope a = dVbg/dVtg = — (etdb)/(ebdt) — —7.72, where eb(t) and db(t) 
refer to the relative bottom(top)-gate dielectric constant and bottom(top)-gate dielectric 
thickness, respectively (Fig. [H^c)). Using ^b = 3.9 (Si02), db = 300 nm and dt = 40 nm , 
we obtained a relative dielectric constant for hBN St — 4.0 which is in good agreement with 



the reported value for bulk hBN 



ll| . The introduction of Vbg and Vtg capacitively changes 



the carrier density in BLG and hence shifts its chemical potential (/i). For a given Vtg, R 
attains a peak value whenever /i is shifted back to CNP by tuning Vbg while giving a finite 

— * — * 

D on BLG that grows in magnitude with Vtg. The total unscreened displacement field D 
on BLG can be calculated using D = [£:b(Vbg — Vbgo)/db — £^t(Vtg — Vtgo)/dt]/2, where Vbgo 
and Vtgo equal 5.84 V and 0.75 V, respectively, determined from the linear fit shown in Fig. 
[U^d). Therefore, the large increase in the resistance peak is a direct consequence of the band 
gap opening (right panel in Fig. W^c)) due to the inversion symmetry breaking by D, which 
appears to be most dramatic in BLG Q. 

At Vtg = and 15 K, R attains a peak value at Vbg — 5.84 V shown as thick red line in 
the upper panel of Fig. [2]^a). For Vtg 7^ 0, the double-peak feature emerges and results from 
the partial coverage of the top-gate region as mentioned earlier. It is then straightforward 
to express R = Re+Rne, where Re(ne) refers to the resistance contribution from top-gate 
(non)effective region. As Vtg increases, the peak value of Re grows rapidly with its position 
moving to higher Vbg value. At Vtg = -8 V, the increase is nearly 9-fold and tends to grow 
further at higher Vtg value. Rne, on the other hand, shows relatively weak dependence on 
the Vtg and can be extracted unambiguously (orange dashed line in Fig. EJ^a)). The sheet 
resistance for the top-gate effective region Rne = ReW/Lc, where W is the width of BLG, 
can then be determined as shown in the upper panel of Fig. Mjo)- At Vtg = 0, the ratio 
Re/R ~ 0.58 is close to the length ratio of Le/Lg ~ 0.59 as expected (Fig. [T](a)). 

The thick red curve in the lower panel of Fig. [2t^a) shows the thermoelectric signal Vg 
as a function of Vbg at Vtg = 0. Vg is nearly zero at CNP and exhibits ambipolar feature 
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where electron-type (Vs <0) and hole-type (Vs >0) carriers can be readily tuned by Vbg. It 
increases rapidly in magnitude with Vbg attaining local extremes and falls down at higher 
Vbg. Similar to the analysis for resistance, we denoted Vs = Vge + Vsne, where Vse(sne) is 
the thermoelectric signal from the top-gate (non)effective region. When turning on Vtg, 
Vse separates from the thick red curve with local extremes occurring at higher Vbg as Vtg 
increase in magnitude as shown in lower panel of Fig. [2]^a). The blue curves are Vg signals 
at Vtg = ±10V that we used for the extraction of the Vsnc shown as the dashed orange line 
in the lower panel of Fig. Mjo)- However, the ratio of Vsc/Vs ~ 0.37 at Vtg = turns out 
to be about 20 % less than the expected ratio of Lc/Lg considering the case of an uniform 
(— VT) across BLG. This implies a smaller temperature gradient (— VT)o under the top-gate 
effective region due to a likely heat shunt through the top-gate(Ti-Au)/hBN. By assuming 
that TEP at zero top-gate voltage is the same throughout the BLG, we can then determine 
ATe and hence TEP Se = Vge/ATg under the top-gate effective region which is shown in 



the lower panel of Figj2]^b) 12[. At 15 K, the difference between local extremes, denoted as 
ASm (lower panel in FigJ^b)), equals 20 /iV/K and is enhanced by more than 4 folds to ~ 
95 /iV/K and 80 /iV/K at Vtg = -7 V and 8 V, respectively. 

The peak values Rneo in log-scale are plotted as a function of the corresponding D values 
at different temperatures shown in Fig. [3](a). The solid symbols that extend to higher D 
are obtained from the same device after cooling down to the base temperature again, which 
shows consistent behavior. Rneo increases exponentially with |D| giving Rneo — 230 kQ at 
D ^ 1.2 V/nm and 15 K, which is nearly a 40-fold increase compared to its value at D ~ 
V/nm. However, we notice that the increase in Rneo drops slightly for D > 0.8 V/nm. 
In order to gain further information on the band gap Eg, we plot the relative conductivity 
Rneo(D = 0)/Rnoo as a function of 100/T for eight different D values such that effects other 
than D can be excluded. The data points are then fitted in a temperature range of 50 
K-300 K using Rneo(D = 0)/Rnco = Aexp(— Eg/2kBT) -|- C, where A and C are constants 
independent of T. The extracted Eg vs. D is shown in Fig. Et^c), which is practically linear 
for D > 0.3 V/nm giving a Eg ^ 100 meV at D ^ -0.9 V/nm close to the value given by 
infrared microspectroscopy [6]. The data points show good agreement with the calculated 
band gap (solid line in Fig. [3]^c)) using a self-consistent tight-binding model. The fitting 
function is based on a simplified model of a narrow band-gap semiconductor with impurity 
band that may originate from disorder and impurity in BLG. The parameter A reflects the 
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D dependence of carrier mobility fic = er/m*. It turns out to increase with |D| (inset of Fig. 
[3]^c)), where the charge re-distribution (screening) and also change in band structure may 
play a role. 

We defined Sm = ASin/2 which is the maximum value of Sg. The increase of Sm relative 
to its value at D = 0, [Sm(D)/Sm(D = 0)] — 1, is shown as open symbols in Fig. |4](a). 
[Sm(D)/Sm(D = 0)] — 1 grows larger at lower temperatures and exhibits a minor asymmetry 
with respect to D. At 20 K, [Sm(D)/SiJT) = 0)] - 1 reaches a value of ~ 4.2 at D = 0.7 
V/nm. According to the Mott relation 16|, TEP can be described by 



^" 3 e cr^9e^^=^' 

where a is the electrical conductivity. Using cr=l/Rne, we can then deduce TEP from the 
sheet resistance data using Sg = ^"7"[^^(|v^)](^^)e=M according to Eq. [H The increase 
of the term [^^(f^^)] with respect to its value at D = is shown as the solid lines in Fig. 
m^a) without any scaling or shifting on them. For T < lOOK and |D| < 0.3 V/nm, we find 
surprisingly good agreement in the solid lines to the measured values of [Sm(D)/Sm(D = 0)]- 
1 (symbols). However, the solid lines starts to deviate downward from the symbols as |D| > 
0.4 V/nm. Based on the Mott relation, the discrepancy should reflect the D dependence 
of the term (dVbg/de)|e=^ that is proportional to the DOS at fi. We argue that, as |D| 
increases, more states near CNP are pushed aside forming a band gap that grows wider 
with |D|. The DOS at the bottom(top) of the conduction(valence) band is then expected to 
increase with increasing |D|. As shown in the inset of FigHl^a), AVbg, which is the width of 
the local extremes for measured Se, equals ~ 12 V at 15 K corresponding to a carrier density 
Um ~ 1 X lO^^cm"^ at which Sm occurs. The gradual increase of at higher |D| is consistent 
with the scenario described earlier that may also partly account for the D dependence of 
A in the inset of Fig. M^c). However, for 200 K and 300 K, the solid lines deviate upward 
instead. The failure of the Mott relation near CNP at high temperature has also been 



191]. It was attributed to the violation for 



reported previously by several authors 
the criteria of keT/ep <^ l(eF is the Fermi energy), which becomes more pronounced in a 
cleaner BLG. We note that the enhancement in TEP is mainly associated with the change 
in band curvature due to D rather than Um that shows relatively weak variation with D. 

The temperature dependence of Sm at different D values is plotted in Fig. |4](b). It 
exhibits a maximum near 100 K giving a value of Sm(D = 0.7V/nm) ~ 180/xV/K. Below 
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50 K, we remark that Sm(D = 0.7 V/nm), which tends to increase even further at higher 
|D|, is comparable to the reported large TEP in sodium cobaltates NaxCo02 (x=0.97 and 
0.88) [2^ and Bi2Te3-based alloys (CsBi4Te6 [21] and Bi2Te3 23|) shown as dotted lines. 
Unfortunately, the detailed information on thermal conductivity n in dual-gated BLG is 
absent for the determination of thermoelectric figure of merit ZT= S^ctT/k, where S, a, 
T and k are thermoelectric power (TEP), electrical conductivity, temperature and thermal 



conductivity, respectively 



22l |. Nevertheless, it was recently pointed out that k in encased 



few-layer graphene can be orders of magnitude smaller due to the quenching of flexual 
phonon mode 
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26l |. which makes dual-gated BLG device a potential candidate for large 



ZT. 

Theoretical calculations including the screening effect are based on Kubo's formula of the 
linear response coefficients jl3|. In the clean limit, the predicted relative increase of S^^i at 
15 K and 300 K are shown as dotted and dashed lines in Fig. ID^a), respectively, which is 
more than three-fold larger compared to our experimental results. The corresponding Um is, 
however, an order of magnitude smaller than the experimental data. We suspect that the 
presence of charge puddles near CNP can be an important factor, which has been revealed 
from scanning tunneling microscopy in exfoliated graphene 27|, l28|] and BLG 291] on Si02 
substrate. The effect of electron- hole puddle gives rise to a large charge inhomogeneity 



5n ~ lO^^cm ^ 
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341 ] near CNP that falls in the same order of magnitude as 



Urn in our device . The coexistence of electron-type and hole-type carriers at close to 



CNP is reminiscent of the finite minimum conductivity 3014331] and also phonon-anomaly 



in bilayer graphene 



34| . This can also cause a significant compensation in TEP near CNP 
as suggested in our experiment. We also remark that the relative increase of Sm depends 
critically on Um as shown in the inset of Fig. Hl^b) obtained from theoretical calculation for 
[Sm(D = 0.6V/nm)/Sni(D = 0)] — 1 at different The value of Um inevitably incorporates 
the impurity doping giving rise to larger impurity scattering at higher Um and hence a smaller 
enhancement in TEP. 

In conclusion, we demonstrate the electric-field tunable band-gap and TEP in dual-gated 
BLG device, which may offer a new platform for innovative science and engineering. Un- 
fortunately, large enhancement in TEP only occurs below 100 K in our devices most likely 
related to the compensation from electron-hole puddles near CNP. It gives Sm(15K, D = 
0.7 V/nm)=48 /iV/K that is comparable to or exceeding existing records for low-T thermo- 
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electric materials. Larger enhancement, in principle, can be realized in a cleaner dual-gated 

— * 

BLG device at higher D. With the advantage of full electric field control on TEP and also its 
carrier polarity, dual-gated BLG device with a proper design can be a promising candidate 
for low-temperature thermoelectric application. 
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FIG. 1: (color online) (a) An optical image of a dual-gated BLG device with shaded area referring 
to the top-gate effective region. The measurement set-up is illustrated with labels in the figure and 
polarity definitions for Vtg and Vbg. (b) A contour plot for the resistance (log-scale) as a function 
of Vtg and Vbg at 200 K. The left panel of (c) illustrates the application of displacement field D 
on BLG using dual-gated geometry (not to scale). The right panel of (c) shows the band gap (Eg) 
opening due to the inversion symmetry breaking by D. (d) The (Vtg, Vbg) values for resistance 
peaks in (b) at 15 K, 50 K, and 200 K. The linear fit shown as red line gives Vbgo = 5.84 V and 
Vtg = 0.75 V. 
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FIG. 2: (color online) The upper panel in (a) shows R vs. Vbg for the dual-gated BLG at 15 K and 
different Vtg values ranging from -10 V to lOV. The lower panel plots the corresponding thermo- 
electric signal Vs at the same temperature. The dashed orange line represents signal contribution 
from the top-gate noneffective region extracted from the blue curves where Vtg = it 10 V. (b) Vbg 
dependence of extracted sheet resistance Rne (upper panel) and TEP Se (lower panel) for top-gate 
effective region. The thick red lines are the signals at Vtg = 0. The definition of ASm is shown in 
the lower panel of (b). 
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FIG. 3: (color online) (a) Rdco (log-scale) vs. D for the dual-gated BLG at 6 different temperatures 
ranging from 15 K to 300 K. Rdco exponentially grows with D. (b) The relative conductivity 
Rnoo(D = 0)/Rnco(D) vs. 100/T at 8 different D values. The red lines are the fitted curves using 
Rnoo(D = 0)/Rneo(D) = Aexp(— Eg/2kBT) + C. The extracted band gap Eg and A parameter as 
a function of corresponding D are shown in (c). The solid line is the calculated band gap using 
self-consistent tight-binding model. 
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FIG. 4: (color online) (a) [Sm(D)/Sm(D = 0)] - 1 vs. D for the dual-gated BLG at 6 different 
temperatures ranging from 15 K to 300 K. The solid lines are the relative increase of [^~( gv"" 
obtained from the Rne data. The dotted and dashed lines are the theoretical prediction at 15 
K and 300 K, respectively. The inset figure plots the width AVbg and its corresponding carrier 
density at which Sm occurs, (b) T dependence of for the dual-gated BLG at 9 different 
D values ranging from -0.8 V/nm to 0.7 V/nm shown as symbols. The TEP values for NaxCo02 
(x=0.97 and x=0.88), Bi2Te3 and CsBi4Te6 are shown as dotted lines for comparison. The inset 
figure shows the theoretical calculation of [Sm(D = 0.6V/nm)/Sm(D = 0)] — 1 vs. n^. 



12 



